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We demonstrate an area law bound on the ground state entanglement entropy of a wide class of 
gapless quantum states of matter using a strategy called local entanglement thermodynamics. The 
bound depends only on thermodynamic data, actually a single exponent, the hyper-scaling violation 
exponent 0. All systems in d spatial dimensions obeying our scaling assumptions and with 0 < d— 1 
obey the area law, while systems with 0 = d — 1 can violate the area law at most logarithmically. 
We also discuss the case of frustration-free Hamiltonians and show that to violate the area law 
more than logarithmically these systems must have an unusually large number of low energy states. 
Finally, we make contact with the recently proposed s-source framework and argue that 9 and s are 
related by s = 2 e . 


When a quantum many-body system is at or near zero 
temperature, quantum entanglement between the con¬ 
stituents leads to qualitatively new phenomena (for a re¬ 
view and references, see e.g. [1, 2]). Although direct mea¬ 
surement of many-body entanglement awaits advances in 
quantum state engineering, it is already an invaluable 
theoretical tool for diagnosing the physics of quantum 
many-body systems. An important role is played by the 
entanglement entropy of a spatial region, which quanti¬ 
fies the amount of bipartite entanglement between the 
region and its complement. 

Here we demonstrate that the ground state entangle¬ 
ment entropy obeys an area law for a wide variety of 
gapless and scale invariant phases of matter. While the 
area law is widely believed to hold for gapped phases 
(e.g. [3]), certain gapless phases are definite exceptions, 
notably conventional metals [4-7] and conformal field 
theories (CFTs) in 1 + 1 dimensions [8]. In addition to 
demonstrating an area law, our results expose a deep 
connection between thermodynamic properties and en¬ 
tanglement properties in scale invariant phases of mat¬ 
ter. Such a connection brings entanglement closer to ex¬ 
perimentally accessible probes. Previous work relating 
the scaling of geometric entanglement to thermodynam¬ 
ics and energy fluctuations includes [9, 10]. 

Given a bi-partite quantum state \?Paa) we study 
the reduced state of the A subsystem, pa = 
aa)aa\) i an d i n particular its entanglement en¬ 
tropy, S(A) = -XY A (pA\og(pA)). When | if) AA ) is the 
ground state of a local Hamiltonian the entanglement 
entropy often obeys an area law, 5'(A) ~ cM meaning 
the entropy is proportional to the size of the boundary of 
A. An area law’s worth of entanglement always appears 
due to short-distance correlations, but to have more than 
an area law’s worth of entanglement intuitively requires 
some long-range structure in the quantum state. For 
example, conventional metals with their Fermi surface 
of low energy electronic excitations violate the area law 
with a multiplicative logarithmic correction [4-7]. The 
intuition for this violation is that metals have a great 
many spatially extended low energy excitations and this 
plethora of low energy modes leads to a plentitude of 


long-range entanglement in the ground state. It is the 
purpose of this paper to make this intuition more precise 
and to extend it to other kinds of scale invariant states 
of matter. 

The argument appeals to ordinary thermodynamic 
properties of the state of matter to quantify the num¬ 
ber of low energy excitations. The connection between 
thermodynamics and entanglement then proceeds by re¬ 
casting the entanglement entropy problem as a problem 
of local thermodynamics, that is, thermodynamics with 
a locally varying temperature. In fact, given a state p A 
arising from a scale invariant ground state, there is an¬ 
other state a A of a local thermodynamic form, whose 
entropy bounds that of p A - 

Assuming the thermal entropy scales with temperature 
like s(T) ~ with 9 the hyperscaling violation ex¬ 

ponent and z the dynamical exponent, we show that the 
corresponding ground state obeys the area law provided 
9 < d — 1 and z is positive and finite. The case 9 = d — 1 
(which occurs for conventional metals) is marginal and 
leads to a multiplicative logarithmic violation to the area 
law. Physically the only assumptions are translation in¬ 
variance, scale invariance, and the validity of local ther¬ 
modynamics within a derivative expansion; the precise 
assumptions are stated just below. 

Problem setup. Consider the ground state | g) = | ipAA) 
of a local d-dimensional Hamiltonian H = H x defined 
on L d sites (energy scale J, range i) with a l/poly(L) 
gap which supports scale invariant physics. Let p A = 
t r ^4 (I <7) (<?|) be the reduced density matrix of region A, 
and let a A denote the maximum entropy state consistent 
with the expectation values of all the H x contained in 
A. In other words, a A is the state of maximum entropy 
which gives the same expectation values as p A for all the 
terms in the Hamiltonian contained in region A. 

Since p A is consistent with its own local data, it follows 
that S(a A ) > S(p A ). Furthermore, a A is a local Gibbs 
state: 

^ = e x P(-Exg a h x / t ( x )) ^ 

where again the H x denote the terms in the local Hamil¬ 
tonian and the T(x ) are constants adjusted so that 
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tr (paH x ) = tr(crAH x ). This local Gibbs form arises from 
maximizing the entropy subject to the constraint that lo¬ 
cal data is correctly reproduced [11, 12]. 

Let Ha denote the restriction of H to region A. Ha is 
a sum of local terms residing within A, Ha = J2 x eA Hx- 
The locality of H implies that ti '(HaPa) is within J\dA\l 
of the ground state energy of Ha (where J = | H x | mea¬ 
sures the coupling strength). The bound arises because 
the lost terms in H are localized at the boundary dA , 
so pA is a ground state of Ha up to excitations near the 
boundary. 

As a local Gibbs state, a a is essentially a thermal state 
of Ha with a position dependent temperature T{x). The 
form of T(x) is tightly constrained since the expecta¬ 
tion value of H x is independent of x in the ground state 
and thus in the state a a as well. Furthermore, the local 
temperature approaches zero away from the boundary to 
ensure that the average excitation energy relative to the 
ground state of Ha is localized near the boundary of A. 
Hence a picture emerges wherein the entropy of a a is 
concentrated near dA. We now make this idea sharp. 

The crucial observation is that we can estimate the 
scaling of the energy and entropy of a a using local (in¬ 
tensive) thermodynamic expressions. In other words, if 
e(T) and s(T) are the bulk thermodynamic energy and 
entropy densities respectively, then for purposes of study¬ 
ing the scaling with region size, we may make the replace¬ 
ments 


H(Ha<Ja) ~ Eg t A + J d d xe(T(x)), (2) 

and 

- tr(<7Alog(oyi)) ~ J d d xs(T(x)), (3) 

where T(x) is the local temperature (inverse of the local 
coefficient of H x in — log(cr J 4 )). As far as bounding the 
entropy is concerned, local thermodynamics is a good ap¬ 
proximation because neglecting correlations between dis¬ 
tant regions (that is, treating the system locally) should 
only increase the effective entropy. Later we give a more 
careful justification of the assumption of local thermody¬ 
namics. 

To use this assumption we must specify the thermo¬ 
dynamic properties of the scale invariant phase. The 
relevant thermodynamic scaling data are the dynamical 
exponent z and the hyper-scaling violation exponent 8. 
z relates energy to momentum as ui ~ k z . 8 controls 
the scaling of the entropy density with length: T~ * is 
a length and s is assumed to scale as d — 8 powers of 
inverse length, so s(T) ~ . 8 controls the extent to 

which the naive scaling with density fails for the entropy 
(and other thermodynamic quantities), hence its name. 
In systems where 8 is non-zero, other fixed microscopic 
length scales make up the units of the entropy density. In 
the example of a Fermi surface, the units are made up by 
8 = d — 1 powers of the Fermi momentum. The scaling 


of the energy density is determined by thermodynamics 
to be e(T) ~ Ts(T). 

An important aspect of the resulting theory of entan¬ 
glement thermodynamics is that it is characterized by 
the scaling exponents z, 8 of the fixed point H itself. 
In particular, the effective temperature T(x) must vary 
smoothly with a; in a way controlled by z and 8 to give a 
translation invariant expectation value for H x (although 
oscillations with other fixed length scales when 0 / 0 
are not ruled out). Remember also that if the subregion 
A grows to encompass the total system then the exact 
ground state is recovered [11, 12], so it must be true 
that T{x) decays as x moves away from the boundary 
of A. For future reference, note that for frustration-free 
Hamiltonians the above analysis further simplifies since 
T(x) = 0 identically. 

Entanglement entropy bound. Now we estimate the 
entropy in a simple geometry. Suppose A is a half-space 
in d dimensions with translation invariance in d — 1 
transverse dimensions. We compactify these transverse 
directions to have size R. Translation invariance implies 
that T(x) depends only on x, the distance from the 
boundary. We also introduce a short distance cut-off, a, 
which could be the lattice spacing, and a long-distance 
cutoff, w, which could be the width of the half¬ 
space (making it actually a long strip of length R^> w). 

// 

T 



The scaling of energy with momentum (inverse length) 
and the absence of any other scale determines the local 
temperature T(x) to be 


T{x) ~ x~ z . 


( 4 ) 


In fact, two other solutions consistent with scale invari¬ 
ance are T = 0 and T = oo, and while T = oo can be 
ruled out (too much energy), T = 0 is actually relevant in 
the context of frustration-free Hamiltonians as discussed 
below. For the special case of Lorentz-invariant systems, 
this relation is rigorous, since the entanglement Hamilto- 
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nian for a half space (= — log Phaif-space) is a boost gener¬ 
ator [13-18]. The scaling form for T(x) is further justified 
in an appendix. Assuming this form for T(x ), the energy 
density is 

e(T(x)) ~ x~ z+e ~ d (5) 

and the entropy density is 

s{T{x)) ~ x e ~ d . (6) 

Integrating the energy density gives 

nW 

ti-{H A a A ) ~ Eg tA + -R d_1 / dxx- z ~ d+e (7) 

«/ a 

which converges as the long-distance cutoff w is taken 
to infinity provided z + d > 1 + 9. For example, any 
relativistic CFT has 0 = 0 and z = 1, so the integral 
converges for all d in this case. Since we must have 
tr(H A a A ) < Eg^ A + J\dA\t we see that any local theory 
which obeys our assumptions must have z + d > 1 + 9. 
Turning to the entropy density, we have 

nW 

-tv(<j A \og(<j A ))~R d - 1 l dxx ~ d+e , (8) 

J a 

so that the integral converges provided d > 1 + 9, inde¬ 
pendent of the value of z\ Since (8) is an upper bound on 
S A , all phases of matter with 9 < d— 1 therefore obey the 
area law. This is our main result. The case 9 = d — 1, 
realized in Fermi liquids, is marginal and gives a loga¬ 
rithmic divergence with w. Since we are upper bounding 
the entropy, it does not follow that phases with 9 = d — 1 
must violate the area law. To violate the area law worse 
than logarithmically, a phase must have 9 > d — 1 (or 
have a large number of ground states). 

Such an inequality was previously deduced [19] for the 
special case of systems with a classical gravity dual. 

Validity of the local approximation. We now give a 
detailed justification for the assumption of local thermo¬ 
dynamics. The main idea is to make the local approxima¬ 
tion better by going to higher “temperature.” A similar 
construction was used for other purposes in [20]. Define 

H a by the equation a A = e Z A and then define 


oa(t) 


e -H A /T 

Z(r) 


( 9 ) 


with Z(t) = tre~ HA ' T . r is a fictitious temperature such 
that t = 1 is the original maximum entropy state a A . 
Ordinary thermodynamics for the local Hamiltonian H A 
implies that the entropy of <r A (r) is a monotonic function 
of r. Furthermore, as r grows, the local approximation 
becomes better and better because the correlation length 
becomes shorter (modulo encountering a classical phase 
transition 1 ). 


1 It would be fascinating to find models where one is forced to 
encounter such a classical phase transition upon raising the tem¬ 
perature, and hence where local thermodynamics need not apply! 


Returning to our entropy estimate above, we now have 
T(x) ~ tx~ z and the local approximation gives 

S(<t a (t)) = S(a A ) ~ _R d_1 / dxx~ d+e . 

J a 

( 10 ) 

Thus for any large but R- and w -independent r, S(t) 
scales with R and w the same way as the local approx¬ 
imation for 5(1) and we expect smaller error from the 
local approximation. 

To give a precise characterization of validity of the local 
thermodynamic approximation, we make use of a deriva¬ 
tive expansion argument familiar from hydrodynamics. 
The local applicability of thermodynamics follows if all 
perturbations vary slowly on the scale of the correlation 
length: if this is so, then the system responds locally to 
the perturbation because the memory of distant regions 
is effectively washed out. In an interacting scale invari¬ 
ant system, the correlation length at finite temperature 
scales as £(T) ~ T~ x l z . Note that free theories, e.g., 
Goldstone bosons [21, 22], may not develop such a corre¬ 
lation length at finite temperature, but being free theo¬ 
ries, we may independently establish the area law (or lack 
thereof). Since the spatially varying temperature is the 
only perturbation, the validity of local thermodynamics 
in the derivative expansion rests on the condition 


*^«i. 


(ii) 


In the half-space geometry with r = 1 we argued that 
T(x) ~ x~ z , so we have £ ~ x and the left hand side 
of the derivative expansion condition is of order one, 
S ~ 1. Hence the derivative expansion is not obviously 
well-controlled in this case. However, we can improve 
the situation by introducing a large but system size inde¬ 
pendent r 1. Then because the derivative expansion 
condition is inhomogeneous in T we have the modified 
condition 



£ vr| 

T 


In the last step, we assumed 0 < z < oo. Thus within a 
formal derivative expansion the large r state should obey 
local thermodynamics. 2 

In an appendix we give two additional results related to 
the validity of the local approximation. The conclusion 
in all cases is the same, but these results give us better 
control over the corrections to the local approximation 
and may be of independent interest. 


2 For an extensive quantity / the derivative expansion implies a 
correction of the form / = fi oca i(l + 0(e~ '^)) with the correc¬ 
tions being due to higher derivative terms which are irrelevant 
and thus do not change the scaling of the entropy. Paranoid 
readers may take 1/5 ~ logf R) to fully suppress the corrections 
at the expense of a weaker bound (unphysical log corrections). 
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Frustration-free Hamiltonians. Here we briefly dis¬ 
cuss the very interesting case of frustration-free gapless 
Hamiltonians. This study is motivated in part by [23] 
which presented an interesting frustration-free Hamilto¬ 
nian defined on a segment of length L with the property 
that the gap above the unique ground state was 
and the entanglement entropy of one half of the segment 
~ log(A). Ref. [24] solves another model that violates 
the area law more than logarithmically, in fact as \fL. 

Consider a frustration-free Hamiltonian H with a 
unique ground state | g). Recall that Ha denotes the 
Hamiltonian restricted to A and pa the state restricted to 
A. By assumption, H is frustration-free, so it can be writ¬ 
ten as a sum of local operators H = )T] X H x where each 
H x is positive and annihilates the ground state. Thus we 
have 

tr ( Pa Ha) = 0 (13) 

since every term in Ha independently annihilates the 
ground state. Hence pa lies within the ground state 
manifold of Ha, and as such, its entropy must be less 
than the maximum entropy state in the ground state, 
the uniform mixture of all ground states. Hence we have 
S(pa) < log(G(if, 4 )) where G(Ha ) is the ground state 
degeneracy of Ha- 

It is interesting to note that the maximum entropy 
state consistent with the local terms in Ha is the (nor¬ 
malized) projector onto the ground state of Ha, 

Pmax = G(Ha) P ha = 0 ' 

Since we may write the ground state projector as 

Ph a =o= Hm e~ pHA , (14) 

( 3^00 

we see that in this case we have a very sharp notion of 
entanglement thermodynamics where the temperature T 
is uniform and zero. 

Now suppose that pa violates the area law. Then it 
must be the case that G(Ha) > e s ^ A \ so Ha has an 
enormous ground state degeneracy which grows with sys¬ 
tem size. Indeed, [23] showed precisely this fact for the 
model with log(L) entropy. Hence we see that these mod¬ 
els must have peculiar thermodynamic properties with 
many low-lying states. This is a precise sense in which 
the frustration-free models deviate from the “reasonable” 
criteria of Ref. [25]. We hasten to add that this in no 
way undermines the interestingness of the frustration- 
free models; we simply understand better now how these 
models differ from the familiar examples of scale invari¬ 
ant states. CFTs on an interval, for example, do not have 
a large ground state degeneracy. In particular, in the 
models of [23, 24], the thermal to entanglement crossover 
function discussed in Ref. [25] must take an unusual form. 


Relation to s-source theory. The above considerations 
show that the thermodynamic exponent 6 plays a crucial 
role in the ground state entanglement properties. This 
connection can be clearly displayed by appealing to the 
notion of s source RG fixed points defined in [26]. The 
idea is to classify ground states based on how much en¬ 
tanglement is necessary to create them. Briefly, a phase 
is an s source fixed point if s copies of the ground state 
at linear size L are necessary to produce one copy of the 
ground state at linear size 2 L. 

All s source fixed points obey the entropy bound 
S(2R) < sS(R ) + fcR d_1 , and we generally expect the 
bound to be saturated. Assuming the bound is saturated, 
we can apply the entropy recursion formula to the half¬ 
space region discussed above assuming that the ground 
state is entangled at all length scales shorter than w (the 
IR cutoff, a correlation length). The entropy is then 


S(R) = k 


R 


d_l log (w/a) 


a ' 


l-l 


E 


(2 d -i)" 


= k 


n—0 

R a ~ 1 (\ — (s/ 2 d ~ i y°s( w / a ) 


a' 


i -1 


1 - (s/2*- 1 ) 


(15) 


Although the coefficients may not be correctly repro¬ 
duced, this expression should correctly predict the scaling 
structure - including the scaling structure of sublead¬ 
ing terms. Taking the logarithm base two and writing 
s = 2 log ( s ) the subleading term scales like 


S . 


sub 


_R d_1 a d -i-iog(s) 
a d ~ 1 ^d-i-logfs) ' 


(16) 


Going back to Eq. 8, the local thermodynamic calcula- 

• • , nd—1 

tion also predicts a subleading term of the form ■ 

Although Eq. 8 applies to the maximum entropy locally 
consistent state a a, the scaling structure of subleading 
terms in the entropy should be identical to those in the 
entropy of the actual subsystem state pa- This claim 
amounts to assuming that there is no phase transition 
encountered in going from pA to a a to (where the 

local approximation is justified). Demanding that the 
subleading terms match gives log(s) = 8 or 

s = 2 8 . (17) 


This is a strong result which establishes an intimate con¬ 
nection between thermodynamic and entanglement prop¬ 
erties. It should be noted, however, that this result may 
not apply to frustration-free Hamiltonians since their lo¬ 
cal entanglement temperature is T = 0. It also remains 
to construct the unitary transformation performing the 
mapping from L to 2 L; we have only shown that the en¬ 
tanglement scaling is consistent with being an s source 


3 We note that [24] describes a modification of these models which 


is not frustration-free but exhibits similar entanglement phenom¬ 
ena; this too must have a large density of low-lying states. 
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fixed point with s = 2®. Examples of explicit RG circuits 
will be provided in [27]. 

Discussion. We gave an argument for the area law 
in a wide class of scale invariant phases of matter. The 
argument worked by mapping to the problem of bound¬ 
ing the entanglement entropy to a problem of estimating 
the entropy of a local thermodynamic state. The key 
piece of thermodynamic data is the hyperscaling viola¬ 
tion exponent 9, while the dynamical exponent z played 
little role. This result firmly establishes the intuition that 
highly entangled states of matter must have many low ly¬ 
ing spatially extended excitations. We also discussed the 
analogous statement for the special class of frustration- 
free Hamiltonians and showed that they must possess 
an anomalously large number of low lying excitations in 
order to violate the area law. Finally, we related our 
considerations to the recently introduced s source theory 
and provided another way to quantify the amount of en¬ 
tanglement in the ground state by relating the parameter 
s to 9. 

The formalism accounts for the entanglement proper¬ 
ties of all experimentally-realized states of matter known 
to us which are stable and reach equilibrium. It also ac¬ 
counts for the entanglement properties of several infinite 
families of candidate states (like conformal field theories 
in general and critical Fermi surfaces) that might describe 
as yet mysterious experimental systems. Requiring sta¬ 
bility seems essential since one can concoct fine-tuned or 
pathological Hamiltonians with highly entangled ground 
states (the frustration-free case likely being an example of 
this). On the other hand, glassy physics and other exam¬ 
ples of failures to come to equilibrium are ruled out not 
because the present formalism necessarily fails (although 
it might), but because those cases deserve a separate and 
careful exposition. Among the states of matter that do 
fall into the present discussion, a number of interesting 
examples are displayed in Table 1. The point of this (not 
exhaustive) table is to emphasize that a great many ex¬ 
perimentally realized quantum states of matter are cov¬ 
ered by the gapped analysis in [26] or the present gapless 
analysis and plausibly fit within the s source theory. 

To elaborate on an interesting case, [29] recently ar¬ 
gued convincingly that the entanglement entropy of a 
diffusive metal (meaning a metal in the presence of static 
disorder with a diffusion pole in the density-density cor¬ 
relator) obeys the area law, unlike its clean cousin, the 
Fermi liquid. This is visible in the present framework as 
follows. The low energy density of states is not strongly 
modified by disorder, so the thermal entropy still scales 
like s ~ T. However, the metal is now diffusive rather 
than ballistic, so energy scales with wavenumber like 
ui ~ k 2 and hence z = 2. Requiring s ~ T^ d ~ e ^ z ~ T 
forces 9 = d — 2 which is below the area law violation 
threshold. 

A more complicated application of the formalism is 
provided by the exciton Bose liquid state of [31]. This 
state of matter defined in d = 2 has a peculiar dispersion 
relation which resembles w 2 ~ k'lk‘i near k = 0. This 


TABLE I: Notation: SSB = spontaneous symmetry break¬ 
ing, I/F QHE = integer/fractional quantum Hall effect, QCP 
= quantum critical point, QED = quantum electrodynamics, 
QCD = quantum chromodynamics. QCD is listed as s = 1* 
because although the theory is ultimately gapped, the correla¬ 
tion length diverges in units of a hypothetical short-distance 
cutoff length which is taken to zero. Some representative 
calculations and further references can be found in [2-8, 20- 
22, 26, 28-30], 


State of matter 

z 

s 

9 

EE 

Insulators, etc. 

Gap 

0 

n/a 

Area 

SSB, discrete 

Gap 

0 

n/a 

Area 

IQHE (invertible) 

Gap 

1 

n/a 

Area 

FQHE 

Gap 

1 

n/a 

Area 

Topological states 

Gap 

1 

n/a 

Area 

SSB, continuous (d > 1) 

1 

1 

0 

Area 

QCP (conformal), d = 1 

1 

1 

0 

Area*Log 

QCP (conformal), d > 1 

1 

1 

0 

Area 

Quadratic band touching 

2 

< 1 

0 

Area 

Fermi liquids 

1 

2 d ~ 1 

d- 1 

Area*Log 

Spinon Fermi surface 

3/2? 

2 d ~ x 

d- 1 

Area*Log 

Diffusive metal, d = 3 

2 

2 d—2 

d — 2 

Area 

QED 

1 

1 

0 

Area 

QCD 

Gap 

1* 

0 

Area 


dispersion relation naively suggests z = 2 and 9 = 0, but 
this is only part of the story. Because the dispersion is 
zero for all k x when k y = 0 (and vice versa) the state 
also possesses lines of zero energy states. The k y = 0 line 
has z = 1 with a variable velocity depending on k x and 
effectively has 9 = 1. Thus although the thermal entropy 
goes like s ~ Tlog(l/T) the entanglement is dominated 
by the 9 = 1 zero energy lines and yields a simple loga¬ 
rithmic violation. This violation has been seen explicitly 
in [32] which realizes the counting argument of [7]. 

Hyperscaling violation necessitates the presence of an 
extra length scale i in the density of states, for example 
in order that the entropy density have units of length _ci : 
s(T) = T^~£~ e . In the case of metals, this length scale 
is the Fermi wavelength. We have assumed in various 
places, in particular in claiming that in the half-plane 
geometry £ = x, and that this microscopic length scale 
does not enter. This assumption is physically sensible 
and accords with all examples we are aware of; a further 
general argument is presented in an appendix. 

The derivative expansion was carefully justified by 
introducing the fictitious temperature parameter r in 
ua(t). This trick gives most directly a bound on the 
entropy of the state of interest, but unless a phase transi¬ 
tion is encountered as a function of r, one concludes that 
the entropy of a a{t = 1) has the same scaling structure 
as the entropy of cta(t 1). Similarly, although the 
entropy of a a is only guaranteed to bound the entropy 
of pa , one again expects the entropies to share the same 
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scaling structure. This is because if A were the entire sys¬ 
tem then a a would equal pA , so we expect that as A is 
made larger the local approximation captures more and 
more of the entropy of pa • This can be explicitly checked 
in some cases, e.g. in Lorentz-invariant systems because 
— log(pA) is local near the boundary of A for any region 
A (and this is where most of the entropy arises according 
to the thermodynamic argument). 

The framework developed here, besides yielding a 
strong argument for the area law in scale invariant quan¬ 
tum states of matter, is of broader interest. Similar ar¬ 
guments have been used successfully for gapped phases 
of matter [26] (there the analog of 9 turns out to be the 
ground state degeneracy). It would be interesting to de¬ 
velop further the idea of entanglement thermodynamics 


into a full fledged theory of entanglement hydrodynam¬ 
ics, e.g. in dynamical settings. For example, does the 
state <7 obey a simple dynamical equation? If so, this 
would be interesting since this equation would naturally 
include dissipation as information is lost into the expo¬ 
nential complexity of the quantum state. 
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Appendix A: Further results on the local 
approximation 

Consider a translation invariant system with thermal 
entropy density s(T). The total entropy of a volume V 
at temperature T is thus Vs(T). Now consider the same 
system with a position-dependent temperature, so that 
the state of the system is 


p{{T{x)}) 


ex P (-Ex#)) 
Z{{T{x)}) 


(Al) 


Suppose now that the temperature T(x) is approximately 
constant, T(x) = T + ST(x) with 5T(x ) <C T. Then the 
state of the system is close to the translation invariant 
thermal state. The first order correction is 

5p = p({T(x)})—p(T) « ^ 1 H x -(H x ))5T(x ^ p(T) 

' ( A2 ) 

We have not been careful about the operator ordering, 
but we will only use this expression within a trace so no 
harm can arise to first order in 5T/T. 

The change in the entropy is 


SS = S(p({T(x)}) - S(p(T)) « -tr (Sp\og(p )). (A3) 


This expression reduces to a two-point function of H x , 
SS « p Y.^ H * H v) - (H,)(H y ))6T(x). (A4) 

x,y 


Translation invariance in the uniform temperature state 
implies that the connected two-point function ( H x H y ) — 
(H x )(H y ) depends only on x — y. Thus after the sum 
over y is performed, the change in the entropy is simply 


where H = J^ x // ' - 

Now compare this expression with the expression from 
local thermodynamics, 

SS « £ »(T(*)) - Vs(T) » •£ d t s(T)5T(x). (A6) 

xGV x 

T8ts(T ) = c(T) is the heat capacity, which is in turn 
related to energy fluctuations by 

c(T) = ^((H 2 ) - (H) 2 ). (A7) 

We conclude that to first order 

s( P anx)})) = Y / mx)) ( as ) 


where s(T ) is the translation invariant entropy density. 
We wish to argue for a stronger statement, namely that 

5(p({T(a:)}))«^ S (T(x)) (A9) 


whenever T(x) varies slowly on the scale of the local cor¬ 
relation length. This is the content of the derivative ex¬ 
pansion discussed in the main paper. To quantify the va¬ 
lidity of the derivative approximation it is useful to look 
at the derivatives of the entropy with respect to T(x). 
It will be useful to distinguish two averages, (...) over 
the full non-uniform state, and (.:)t(x) over the uniform 
state with uniform temperature T(x). 

Computing /i = d T ( x) S(p{{T(x)})) gives 

_/ H x ^ (H x ) v H y \ _ v (H X H V ) - (H x ){H y ) 
71 \ T{xY # T (y)] # T{xYT{y) 

(A10) 

Similarly, computing f 2 = d T(x) E y s{T(y)) gives 


r _ C (T(x)) _ (H x ) T ( x ')(Hy) T ( x j 

( A11 ) 

Comparing these two expressions gives a quantitative 
condition. The connected two-point function of the 
Hamiltonian “density” must be approximately equal in 
the two states. Assuming this connected correlation func¬ 
tion is exponentially decaying (due to the finite correla¬ 
tion length) and also exponentially insensitive to distant 
properties of the system (again, the correlation length), 
one concludes that the entropy of the temperature de¬ 
pendent state is well approximated by the sum of the 
translation invariant entropy density evaluated at the lo¬ 
cal temperature. 


Appendix B: Further justification of effective 
temperature scaling 


6S = ^((H 2 ) - (H) 2 )Y,6T(x) 

X 


(AS) 


On general scaling grounds the effective temperature in 
a a was assumed to go like T{x) ~ l/x z . The only other 











possibility is T = 0 which is realized in frustration-free 
systems and which we do not consider in the appendix. 
However, this argument is slightly subtle when 6 ^ 0 
since dimensional analysis implies that at least one addi¬ 
tional length, call it £, does not decouple in the low en¬ 
ergy limit. Perhaps the effective temperature only scales 
like T(x) ~ f(x/t)/x z with / an undetermined scaling 
function? 

The first argument in favor of / = 1 begins with the 
observation that the (momentum) dimension of the en¬ 
ergy density is d — 9 + z. In the presence of an infrared 
length scale x (like the distance to a boundary), the en¬ 
ergy density should to be modified by an additive shift 
of the form 5e/x d ~ 8+z . This is exactly what is obtained 
above when T(x) ~ l/x z . 

The second argument doesn’t directly constrain / but 
says that either the entropy bound is good or translation 
invariance of the local data must be broken. If / = 1 
then the analysis above goes through and the entropy is 
bounded by an area law for 9 < d — 1. Suppose now 
that / goes to zero as x —> oo. Then a quick calculation 
confirms that the derivative expansion will not be valid 
as the system is approaching the ground state too rapidly 
as x —> oo. However, the effective temperature is much 
lower than in the case / = 1, so the entropy of a state 
with decaying / will be upper bounded by the entropy 
of a state with / = 1 so that the entropy bound is still 
valid. Finally, suppose f goes to infinity as x —► oo. 
Then the derivative expansion will be arbitrarily good 
as x increases, but then it will be possible to detect a 
failure of translation invariance in the local data. This is 
so because a correlation length much smaller than x will 
imply that the boundary at x = 0 is not observable and 
thus the expectation value of H x will depend on x. 

The third and final argument proceeds by analyzing 
energy fluctuations and directly shows / = 1 given the 
assumption that energy fluctuations in the maximum en¬ 
tropy state scale the same way as in the true state. This 
assumption is again made plausible by the fact that the 
true ground state is recovered when A is the total sys¬ 
tem. The argument begins by observing that there are 
two ways to compute energy fluctuations, by directly in¬ 
tegrating the connected energy-energy correlation func¬ 
tion and by integrating the heat capacity over all space. 


Let A h denote the spatial scaling dimension of the en¬ 
ergy density (this dimension can be different from d—d+z 
which is better understood as a temporal scaling dimen¬ 
sion). Then the connected 2-point function (H x H y ) will 
decay like \x — y\ 2AH and the energy fluctuations (fluctu¬ 
ations of Ha) of a strip will have a UV finite term going 
like 


A H\ 


L d ~ 1 

Aff d 1 


(Bl) 


A h may be determined from ordinary thermodynamics 
by observing that energy fluctuations are given by T 2 c(T) 
where c(T) ~ T^ d ~ e ^ z is the heat capacity. At temper¬ 
ature T the energy fluctuations are obtained from the 
k = 0 limit of the Fourier transform of the equal time 
energy-energy 2-point function. Scaling again determines 
that ( HH)(k — > 0) ~ j'(2A H -d)/z w h ere the — d is from 
the integral over space defining the Fourier transform. 
Setting the two expressions for the energy fluctuations 
equal gives 2A H = 2d + 2z — 9. Note a hidden assump¬ 
tion here, that all directions in space scale the same way; 
this assumption can be violated, e.g. in a 2d array of 
decoupled Id wires, but a more general analysis can be 
performed in such cases leading to the same overall con¬ 
clusion. 

Turning now to the local thermodynamic calculation 
for the strip in the ground state, one sees that if and 
only if T ~ l/x z will the energy fluctuations obtained 
by integrating T(x) 2 c(T(x)) be equal to those obtained 
from the double integral of (H x H y ) over A. Hence requir¬ 
ing that the scaling of energy fluctuations be reproduced 
forces T(x) ~ l/x z . 

As a final note, T(x) may also contain terms which 
oscillate at a wavelength set by i. Since t is a micro¬ 
scopic length these oscillations are far too rapid for the 
local approximation to be applicable on that scale. This 
is in fact why such oscillations are allowed (otherwise the 
local approximation would lead to badly non-translation 
invariant local data). If such oscillating terms are sub¬ 
leading compared to the dominant scaling of T(x) then 
all is fine. An envelope to the oscillation decaying slower 
than l/x z is again ruled out because it would lead to 
non-translation invariant local data. 



